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The last unknown contribution to hydrogen energy levels 
at order ma^, due to the slope of the Dirac form factor at 
three loops, is evaluated in a closed analytical form. The 
resulting shift of the hydrogen nS energy level is found to be 
3.016/n^ kHz. Using the QED calculations of the IS" Lamb 
shift, we extract a precise value of the proton charge radius 
r„ = 0.883 ± 0.014 fm. 



Precision experiments with hydrogen and, more gener- 
ally, with hydrogen-like atoms serve as an excellent labo- 
ratory to test theoretical approaches to bound state QED 
(for a recent review see, e.g. Q). These experiments ad- 
dress a number of features of the simplest atoms, such as 
the energy levels of the ground and excited states, and 
the corresponding lifetimes. 

In recent years we have seen remarkable progress in 
the experimental study of the hydrogen atom. In par- 
ticular the accuracy of the IS Lamb shift measurements 
has increased dramatically over the years . All mea- 
surements arc consistent with each other and the most 
accurate value so far was determined in Ref. M : 



AE{lS%^p = 8 172 837(22) kHz. 



(1) 



This result was obtained by analysing the most precise 
measurements for the transition frequencies in hydrogen 
(see for details). 

Theoretically, since the ratio of the electron mass m 
to the proton mass M is very small, it is convenient to 
write hydrogen energy levels as a double expansion in a 
and m/M. The corrections which survive in the limit 
M — > c» are known as non-recoil corrections, the other 
corrections are called recoil ones. It is further convenient 
to organize non-recoil corrections in powers of a"(ZQ!)', 
assigning an auxiliary notation Z for the proton charge. 
In this case the correction a'^iZa)'- describes a contri- 
bution of all diagrams with I — 3 Coulomb photons ex- 
changed between the electron and the proton and with 



* e-mail: 
^e-mail: 



melnikovQslac .stanford.edu 
timo@particle.physik.uni-karlsruhe.de 



n photons emitted and absorbed by the electron. The 
general expression for the 72S-lcvel shift can be written 
as: 



AE — A_Enon-rccoil + AE, 
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where we have also added a contribution of muons and 
hadrons to photon vacuum polarization and a contribu- 
tion due to the proton structure (see a discussion below). 

We begin with non- recoil corrections. These correc- 
tions can be parameterized as (see e.g. ||^-[lO|]): 
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where to^ = mM / (m + A/) is the reduced mass of the 
electron and L = logm/ (mr{Za)-) and the function 
A(Za) contains all higher order terms in the expansion in 
Za. We have also indicated that the higher order correc- 
tions contain a logarithm of the fine structure constant 
in the third power . 

All the terms in the above equation, with the exception 
for C40, are currently known. It is the purpose of this 
paper to report on the calculation of the last missing 
ingredient in C40, the slope of the Dirac form factor at 
zero momentum transfer. 

The hydrogen atom is formed because of a Coulomb 
interaction of the proton with the electron. The interac- 
tion of the virtual photon with the electron on its mass 
shell can be parameterized by the so-called Dirac and 
Pauli form factors: 

w(p2)r^u(pi) u{p2) (i^i(g^)7M + *'^M"^^^2('7^)) ^(Pl)' 

where the u{p) are the electron spinors in the initial and 
final state and q is the momenta carried away by the 
photon. The momenta satisfy the relation q = P2 — Pi- 
An important consequence of QED gauge invariance and 
the electron charge definition is that the Dirac form factor 
equals unity at zero momentum transfer i^i(O) = 1, to all 
orders in the coupling constant. The Pauli form factor at 
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zero momentum transfer describes an interaction of the 
electron spin with the homogeneous magnetic field; it is 
the electron anomalous magnetic moment. 

Let us now turn to the contribution of the Dirac form 
factor to hydrogen energy levels. The typical momenta of 
the Coulomb photon in hydrogen is given by the inverse 
Bohr radius gtyp ~ ma, where m is the electron mass. 
Compared to the electron mass, this momentum transfer 
is quite small and for this reason one can Taylor expand 
the Dirac form factor in powers of q'^ /m?. One obtains: 



where, as we already mentioned, the first term equals 
unity to all orders in the coupling constant. The slope of 
the Dirac form factor, Fl in the above formula, can be 
written as a series in a: 



a check of the calculation, and in particular on our solu- 
tion of the system of recurrence relations, we have repro- 
duced also the analytical value of the three-loop electron 
anomalous magnetic moment. Details of our calculation, 
including the results for the individual diagrams, will be 
presented in a separate publication. 

Our result for the slope of the Dirac form factor is |^: 
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In this paper we consider the first three terms in this 
series. 

Two things should be noted at this point. The first 
and second order corrections to the slope of the Dirac 
form factor contribute to the coefficients A^q and B4Q in 
Eq.(||), respectively. Moreover, the one loop slope ^gi^opg 
is not infrared finite, as can be seen from the fact that 
the term A41 log(a) appears in Eq.(||). This infrared di- 
vergence gets removed if one takes into account that the 
electron is not on its mass shell in the bound state. In the 
language of the effective theories, the one loop slope in 
Eq. corresponds to "hard" contributions to energy lev- 
els. Since the two and three loop slopes of the Dirac form 
factor are infrared finite, the off-shcUness of the electron 
is irrelevant for these contributions. 

Being divergent, the one loop slope depends upon the 
chosen regularization. In the present calculation we have 
used dimensional regularization (the space-time dimen- 
sion is Z) = 4 — 2e) for both ultraviolet and infrared 
divergences. 

Let us briefly describe how the actual calculation of 
^sioDc been done. After applying a projection opera- 
tor for the Dirac form factor on the electron-photon ver- 
tex and after performing the Taylor expansion in the pho- 
ton momentum transfer, one obtains diagrams of the self 
energy type. There are four basic topologies which ap- 
pear in this calculation, characteristic example diagrams 
are depicted in Fig.l. For each of the topologies one 
writes down a system of recurrence relations obtained by 
the use of integration- by-parts |l^. Solving this system, 
it is possible to show that any integral which belongs to 
the above topologies can be expressed through 17 master 
integrals. Luckily, these seventeen integrals have already 
been computed in the course of the analytical calcula- 
tion of the electron anomalous magnetic moment and 
hence can be taken from there. Let us mention, that as 
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FIG. 1. Examples of electron-photon vertex diagrams that 
correspond to the four different integration topologies. 

We now turn to the contribution of the slope of the 
Dirac form factor to hydrogen energy levels. It is easy to 
see that the slope of the Dirac form factor gives rise to a 
perturbation of the Coulomb potential: 

^Kiopc(r) = ^5'{r)Fi. 
This perturbation delivers a nS-level energy shift: 



AEsiope = {i'{r)\5Vsiope{rM{r)) 

Here 1-0(0)1^ = {rriraZ)^ / {ttji^) is the square of the hy- 
drogen wave function at the origin. The correction in- 
duced by the three-loop slope of the Dirac form factor is 
then0: 



^The first correct numerical result on the two loop slope of 
the Dirac form factor was published in Ref. |l4|. The analyt- 
ical result was obtained in jTHj. 

■^In what follows numerical results are given in frequency 
units using AE /\E/{2ttK). 
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where cRoa is the Rydberg constant in MHz. 

Using the values for the Rydberg and the fine structure 
constant lll6| 



cRoo = 3 289 841 960.367(25) MHz, 
a = 1/137.035 999 76(50), 



we arrive at: 



\n) = — 7T- kHz. 
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The contribution due to the slope of the Dirac form fac- 
tor was the last unknown contribution to the hydrogen 
energy levels at order a^{Za)*. The two other contri- 
butions to the coefficient C40 come from the three-loop 
electron anomalous magnetic moment and the three-loop 
vacuum polarization correction to the Coulomb propa- 
gator. These contributions can be extracted from the 
literature [[3|jl|]. 

Taking all three contributions into account, we obtain 
the following expression for the coefficient C40 for the 
S-levels: 
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which results in a three-loop correction to the nS-level 
Lamb shift: 
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In the above equation we have displayed the contribu- 
tions due to the three-loop slope of the Dirac form fac- 
tor, the three-loop anomalous magnetic moment of the 
electron and the three-loop photon vacuum polarization 
function separately to emphasize a strong cancellation 
which occurs between vertex and vacuum polarization 
corrections. Thanks to this cancellation, the correction 
turns out to be quite small numerically. 

We now turn to other results available in the literature. 



For the coefficients Aij and Bij in Eq. 



we user 



^See Ref. [Q] for reference to earlier work. Some of the num- 
bers in Ref. H have been updated according to recent cal- 



culations in Refs. |lq-|2C|] . The value of the non-perturbative 
function A{lS){Za) for Z = 1 is extracted from Ref. for 
the self energy correction and from Ref. for the vacuum 
polarization. 
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The only known non-recoil correction at order 0(ma^) 
is the triple logarithmic enhanced contribution [|ll[ : 
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This correction is included in the central value for the IS 
Lamb shift quoted below. 

Consider now the recoil corrections. Part of these are 
already included in Eq. ^ through its dependence on the 
reduced mass of the electron. The terms which go beyond 
this approximation are: 
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computed in Refs. ||2S 

The vacuum polarization effects due to muons and 
hadrons can be extracted from Ref. ||2^ and give AE'vp = 
-8.5 kHz for the IS Lamb shift. 

The last contribution in Eq.(||) which has been not 
considered so far is the one due to proton structure. The 
major part of this correction to the Lamb shift is param- 
eterized through the proton charge radius but there is 
also a proton self energy correction which goes beyond 
that approximation [p4| . The complete correction reads: 
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There are two different values for the proton charge ra- 
dius 0.805(11) fm and = 0.862(12) fm obtained 
in p6|]27| ], respectively, by analysing the electron proton 
scattering data. The smaller value, — 0.805(11) fm, 
seems to give a serious disagreement between precision 
atomic measurements and the theoretical predictions 
lE9i. For this reason wc do not consider it here. The 
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data of Ref. on the other hand, was recently reanal- 
ized in |2^ , where a normahzation of the proton charge 
form factor at zero momenta transfer was also treated as 
a free parameter in the fit. The analysis of Ref. |^ leads 
to a larger proton charge radius rp — 0.877(24) fm. The 
final value for the IS-shift strongly depends on the value 
of the proton charge radius: 

A£;(lS')thoory = 8 172 778(16) (32) kHz, (14) 
Ai;(lS')thoory = 8 172 819(16)(66) kHz, (15) 

where the values are given for rp — 0.862(12) fm and 
rp = 0.877(24) fm, respectively. The first error in all 
the above equations is the theoretical uncertainty due 
to still uncalculated higher order corrections to energy 
levels beyond the a'' order (see a discussion below). The 
second error is due to the uncertainty in the experimental 
values of the proton charge radius. The uncertainties in 
Rydberg and the fine structure constants are not relevant 
at the present level of precision. 

Comparing these numbers with the most recent mea- 
surement of the IS-level Lamb shift 1711 
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172 837(22) kHz. 
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we conclude that the larger values for the proton charge 
radius seem to give an agreement between the theory and 
experiment. 

An important part of the error in A_Z?(lS')thcory is 
due to uncalculated higher order corrections. For the 
IS' Lamb shift the uncertainty [29| ] was estimated to be 
about 40 kHz. This number is the linear sum of a 8 kHz 
uncertainty in the self energy correction, the 16 kHz un- 
certainty caused by the unknown a^(Za)^ terms and the 
16 kHz uncertainty due to unknown a^{Za)'^ terms. The 
result of Ref. |^l[] removes the first uncertainty and our 
calculation of the three loop slope of the Dirac form fac- 
tor removes the last one. The total uncertainty in the 
theoretical predictions in the IS Lamb shift is therefore 
reduced to 16 kHz. 

Turning the problem around, we note that the small 
theoretical uncertainty on the 15 Lamb shift permits an 
extraction of the proton charge radius by comparing ex- 
perimental and theoretical results. We then arrive at the 
precise value of the proton radius rp — 0.883 ± 0.014 fm. 

In conclusion, we have computed the three-loop slope 
of the Dirac form factor. Thanks to this calculation the 
theoretical uncertainty in the predictions for the IS Lamb 
shift is reduced. Comparison of the theoretical and ex- 
perimental results for the IS* level shift permits an accu- 
rate determination of the proton charge radius. Further 
improvements in theoretical predictions for the IS level 
shift would be possible if subleading a^{Za)^ log^ a cor- 
rections are calculated. Only then can the theoretical 
uncertainty be brought down to several kHz and can the 
potential of the recent measurement of the 18-25 transi- 
tion frequency ^ be fully exploited. 
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